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ABSTRACT 

Elasticity  of  the  liner  and  its  effects  on  cavitation  were  numerically 
determined  for  a  squeeze  film  damper  subjected  to  dynamic  loading.  The 
loading  was  manifested  as  a  prescribed  motion  of  the  rotor  undergoing 
noncentered  circular  whirl.  The  boundary  conditions  were  implemented  using 
Elrod's  algorithm  which  conserves  lineal  mass  flux  through  the  moving 
cavitation  bubble  as  well  as  the  oil  film  region  of  the  damper.  Computational 
movies  were  used  to  analyze  the  rapidly  changing  pressures  and  vapor  bubble 
dynamics  throughout  the  dynamic  cycle  for  various  flexibilities  in  the  damper 
liner.  The  effects  of  liner  elasticity  on  cavitation  were  only  noticeable  for 
the  intermediate  and  high  values  of  viscosity  used  in  this  study. 

INTRODUCTION 

Earlier  work  by  Jacobson  and  Hamrock  (Refs.  1  and  2)  using  high-speed 
photography  revealed  the  existence  of  vapor  cavitation  in  dynamically  loaded 
journal  bearings.  A  polymethylmethacrylate  (PMMA)  housing  was  used  for 
visualization  of  the  cavitation.  In  their  studies  they  investigated  the 
effect  of  roller  material,  the  influence  of  L/D,  effect  of  static  and  dynamic 
eccentricities  and  influence  of  speed.  A  theoretical  model  implementing  the 
Elrod  algorithm  (Ref.  3)  was  developed  by  Brewe  (Ref.  4)  for  the  dynamically 
loaded  journal  bearing.  The  algorithm  conserves  mass  throughout  the 
computational  domain  and  essentially  mimics  the  Jakobsson-Floberg  and  Olsson 
boundary  conditions  (Refs.  5  and  6).  Later  some  of  the  earlier  experiments 


were  used  to  corroborate  the  theoretical  model  (Ref.  7).  A  prediction  of 
32  ms  bubble  life  for  nonstationary  cavitation  was  obtained  from  the 
theoretical  model.  Using  high-speed  photography  for  the  same  set  of 
conditions,  a  bubble  life  of  23  ms  was  observed  in  the  laboratory.  This 
discrepancy  was  attributed  to  deformation  and  temperature  effects  as  well  as 
some  uncertainty  in  correlating  actual  life  with  duration  of  visibility.  To 
achieve  a  more  representative  comparison  it  is  necessary  to  include  the 
thermal  effects  in  the  model  as  well  as  deformation  of  the  shaft  and  housing. 
This  analysis  represents  a  first  step  toward  including  first-order  deformation 
effects  with  the  Elrod  algorithm.  As  such  the  model  assumes  that  the  bearing 
is  essentially  a  thin  tube  surrounded  by  a  rigid  housing.  Since  the  periphery 
of  the  bearing  is  much  larger  than  its  thickness,  the  radial  deformations  of 
the  latter  at  a  point  may  be  assumed  to  be  proportional  to  the  pressure  at  the 
point.  This  model  is  commonly  referred  to  as  the  Winkler  foundation  and  has 
been  used  by  numerous  investigators,  i . e . ,  (Refs.  8  to  12).  At  approximately 
the  same  time  that  Dowson  and  Taylor  (Ref.  9)  published  their  work  for  a 
compliant  thrust  bearing  application,  Castelli,  et  a! .  (Ref.  13)  published 
their  findings  using  an  exact  analytical  approach.  In  this  paper  and  in  a 
discussion  to  the  Dowson  paper  (Ref.  9),  the  limitations  of  the  Winkler 
foundation  are  discussed.  Based  on  this,  the  thin  liner  materials  used  in 
this  investigation  appear  to  be  suited  to  the  assumptions  of  the  theory. 
Realistic  edge  effects  are  not  expected  however. 

The  purpose  of  this  paper  is  to  study  the  elastic  effects  of  the  bearing 
liner  on  cavitation  in  a  squeeze  film  damper.  It  has  been  shown  for  steadily 
loaded  bearings  that  increasing  the  liner  flexibility  increases  the  cavitation 
angle  (Refs.  14  to  17)  using  Reynolds  boundary  conditions.  It  is  believed 
that  under  dynamic  conditions  the  elasticity  of  the  iiner  may  influence  the 
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bubble  dynamics  and  hence  the  damping  characteristics  of  the  squeeze  film 
damper. 
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13  liquid  bulk  modulus,  N  m_2 

Y  angular  position  of  minimum  film,  rad 

A  dimensionless  radial  deflection,  S/AR 

AR  radial  clearance,  m 

S  radial  deflection,  m 

e  eccentricity  ratio,  e/AR 

9  fractional  film  content  in  cavitation  zone;  density  ratio  (p/pc)  in  full 

f i lm  zone 

p  dynamic  viscosity,  N  s  m_2 

p  fluid  density  within  full  film  zone,  kg  nr  3 

pc  fluid  density  within  cavitated  zone,  kg  m~3 

cp  angular  coordinate  relative  to  minimum  film  line,  rad 

wq  orbital  angular  velocity  of  journal  center  about  fixed  point  relative  to 

housing  center,  rad  s-1 

us  angular  velocity  of  journal  about  its  own  center,  rad  s'1 

THEORY 

Cavitation  Algorithm 

A  computer  program  was  developed  based  on  Elrod's  algorithm  (Ref.  3)  that 
combines  a  control  volume  approach  to  deriving  the  finite-difference  equation 
with  a  switching  function  that  automatically  handles  the  moving  boundary 
condition  at  the  liquid-vapor  interface.  The  complication  in  satisfying  the 
boundary  condition  is  avoided  by  defining  a  single  independent  parame^r  9 
which  links  the  full-film  region  and  the  cavitated  region.  This  parameter  is 
defined  by  the  ratio  of  the  local  density  to  that  of  the  cavitation  density, 
i  .e. , 
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It  follows  that  in  the  cavitation  region  0  <  1.0  and  In  the  full-film 
region  0  >  1.0.  A  switch  function  g  was  defined  as  "on"  (g  =  1)  and  "off" 

(g  =  0)  in  the  full-film  and  cavitation  regions  respectively.  By  using  this 
property,  a  universal  Reynolds  equation  was  defined  as  indicated: 

^  +  $  •  m  =  0  (1) 

where  m  represents  the  lineal  mass  flux  and  is  given  by 

—  =  ^  |3g(0)$r  (2) 

Pc  2  1 2y 

In  the  cavitation  region,  the  second  term  vanishes  since  g  =  0.  This 
equation  was  integrated  in  time  by  using  an  Euler's  method  to  advance  in 
time.  At  every  time  step  the  difference  equation  was  solved  by  an  Alternating 
Direction  Implicit  (ADI)  scheme.  The  switch  function  was  updated  after  each 
half-time  step.  Occasionally,  during  the  collapse  of  the  vapor  bubble,  liquid 
was  piling  up  at  the  boundary.  This  was  due  to  an  inherent  time  lag  of  the 
switch  function  with  the  calculation.  Consequently,  the  switch  function  was 
immediately  updated  to  accommodate  the  movement  of  the  boundary,  and  the 
calculation  was  reiterated. 

The  results  of  the  program  developed  were  authenticated  by  comparison  to 
an  experiment  of  Ref.  6.  When  compared  to  the  pseudo-Gumbel  boundary 
conditions  (Ref.  4),  the  results  predicted  a  significant  difference  in  the 
computed  load  capacity  (up  to  20  percent).  Furthermore,  it  was  found  that  the 
size,  life,  and  location  of  the  cavitation  bubble  were  considerably  affected 
by  imposing  the  pseudo-Gumbel  boundary  conditions. 

Fluid  Velocity  Under  Dynamic  Loading 

Consider  the  bearing  geometry  shown  in  Fig.  1.  As  shown,  the  minimum-film 
thickness  rotates  clockwise  around  the  bearing  at  a  whirl  frequency  of  y. 
Jacobson  and  Hamrock  (Ref.  2)  show  that  from  this  diagram  y  can  be  evaluated 
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Y  =  tan' 


ed  sin(u)dt) 

e„  +  c  .  COS(co.t) 
s  d  d 


so  that 


e  ,[e„  COS(oi.t)  +  c  .]  a), 

d  s  d  d  d 

Y  =  - 2 - 

e„  +  e  ,  +  2e„e  .  COS(oj.t) 
S  d  S  d  d 


The  sliding  velocity  of  the  journal  is 

U  =  (0J5  -  2  y)R 
Film  Thickness  and  Deformation 

The  film  thickness  for  a  bearing  having  a  rigid  liner  can  be  given  as 

h  =  AR  -  e  cos  <p 

where  the  total  eccentricity  e  defined  as  the  instantaneous  distance 


between  Ob  and  Oj  is 


ed  sin  (todt) 
sin  y 


If  one  considers  deformation  of  the  bearing  surface,  the  film  thickness 


can  be  written  as : 


h  =  AR  -  e  cos  <p  +  & 


or  normalizing  with  respect  to  the  clearance,  AR,  we  have 

H  =  1  -  e  cos  cp  +  A  (5 

where  6  represents  the  radial  deflection  and  A  =  S/AR. 

This  study  considers  only  1st  order  effects  of  the  elastic  deflection  of 
a  soft  bearing  1 >ner  on  damper  performance.  These  deflections  are  calculated 
using  the  assumptions  adopted  by  Higginson  (Ref.  8)  for  a  soft  liner  material 


rigidly  enclosed,  i .e . , 


A  =  Bp 


This  simple  expression  can  be  derived  by  assuming  the  liner  to  be  an 
elastic  material  that  deforms  linearly  due  to  local  pressures  under  conditions 
of  plane  strain.  Accordingly,  the  parameter  B  can  be  shown  to  be 
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Since  the  film  thickness  is  directly  proportional  to  the  local  pressure, 
the  solution  of  Eq.  (1)  requires  nonlinear  techniques.  An  iteration  procedure 
was  introduced  that  required  the  sum  of  the  pressure  field  to  deviate  by  no 
more  than  0.05  percent  due  to  iterations  in  film  thickness. 

GEOMETRY  AND  OPERATING  CONDITIONS 

The  damper  (bearing)  simulated  in  this  study  is  based  on  an  experimental 
study  of  a  journal -bear i ng  squeeze  film  damper  submerged  in  Dexron  lubricant 
(Ref.  7).  The  range  of  operating  conditions  are  listed  in  Table  I. 

RESULTS  AND  DISCUSSION 
Effect  of  Flexibility  on  Cavitation 

The  effect  of  flexibility  on  vapour  cavitation  was  studied  for  a  journal 
undergoing  noncentered  circular  whirl.  Three  liner  materials  were  used  in 
this  investigation  to  achieve  the  range  of  flexibilities.  These  materials 
were  steel,  polymethylmethacrylate  (PMMA)  and  polytetraf 1 uorethyl ene  ( PTFE ) . 
Table  II  lists  these  materials  and  their  elastic  constants  in  order  of 
increasing  flexibility  for  a  range  of  different  viscosities.  The  percent  of 
time  that  a  vapor  bubble  is  predicted  for  one  full  orbit  of  the  journal  is 
listed  under  "Bubble  life."  To  the  right  of  this  column  are  listed  bubble 
inception  and  collapse.  These  values  are  given  in  milliseconds  of  time 
relative  to  the  moment  that  the  minimum  eccentricity  of  the  orbit  occurs. 

Thus  at  startup,  e  =  0.1  and  t  =  0.0.  The  orbital  period  of  the  journal  is 
67.7  ms  during  which  time  the  eccentricity  ratio  varies  between  0.1  and  0.9. 
Note  that  as  the  viscosity  increases  the  bubble  life  increases  to  the  extent 
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that  cavitation  exists  throughout  the  entire  cycle  for  both  the  steel  liner 
and  the  PTFE  liner.  However,  the  effect  of  increased  flexibility  retards  the 
sensitivity  of  bubble  life  to  increases  in  viscosity.  This  appears  to  be 
attributed  to  the  mechanism  leading  to  the  collapse  of  the  bubble.  From  the 
table  we  see  that  for  the  steel  liner  the  bubble  collapse  is  delayed  as 
viscosity  is  increased.  But  for  the  PTFE  liner  the  bubble  collapse  appears 
insensitive  to  viscosity  changes. 

From  the  standpoint  of  a  viscous  damper,  it  is  logical  that  cavitation 
can  only  serve  to  decrease  the  damping  characteristics.  The  above  results 
lead  one  to  infer  that  as  the  viscosity  is  increased  better  damping 
characteristics  can  be  achieved  using  a  PTFE  liner  over  a  steel  liner  for  a 
limited  range  of  viscosities.  Further,  increasing  the  viscosity  should  give 
better  damping  characteristics,  but  the  accompanying  increased  bubble  life 
could  deter  that  desirable  effect. 

Cavitation  and  Pressure  Distribution 

Figures  2  and  3  illustrate  the  responding  pressure  variation  to  the 
motion  of  the  shaft  undergoing  noncentered  circular  whirl  for  two  contrasting 
situations.  Figure  2  represents  the  events  taking  place  for  a  rigid  (steel) 
liner  with  parameter  B  =  0.000038.  Figure  3  contrasts  those  same  events  for 
a  soft  (PTFE)  liner  having  B  =  0.004440.  In  both  cases,  the  viscosity  was  a 
constant  0.099  Pa  •  s.  In  each  figure,  part  (a)  illustrates  the  starting 
position  (i.e.,  e  =  0.1).  Parts  (a)  to  (f)  depict  several  key  events  in  the 
clockwise  circular  motion  of  the  journal  center  within  the  bearing. 

Figure  2(a)  represents  the  position  of  the  journal  within  the  housing  and 
the  associated  pressure  distribution  at  the  instant  that  the  eccentricity  is 
at  a  minimum.  The  dark  oblong  region  indicates  the  shape  and  position  of  a 
vapor  bubble  downstream  of  the  minimum-film  line.  The  position  of  the  bubble 
within  the  small  clearance  space  of  the  journal  bearing  is  also  indicated.  At 
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this  particular  instant  the  tendency  of  the  pressure  flow  along  the  side  of 
the  bubble  to  cause  the  bubble  to  collapse  becomes  greater  since  it  is 
proportional  to  h^.  The  journal  center  moves  in  a  clockwise  orbit  from 
Figs.  2(a)  to  (b)  (nearly  one-tenth  of  a  revolution).  While  the  eccentricity 
is  increasing,  the  minimum  film  is  decreasing  which  gives  rise  to  the  squeeze 
action  in  the  vicinity  of  the  minimum-film  line.  A  slight  rise  in  the 
pressure  field  can  be  observed  there.  Note  in  Fig.  2(b)  the  collapse  of  the 
bubble  occurs  ~6  ms  later  than  in  3(b). 

Cavitation  is  absent  from  Fig.  2(b)  to  (c).  In  Fig.  2(c),  a  new 
cavitation  bubble  appears,  and  the  pressure  generation  in  the  region  of  the 
minimum-film  line  becomes  significantly  large  because  of  the  squeeze  effect 
taking  place  there.  Figure  2(d)  represents  the  configuration  of  the  bearing 
when  the  maximum  pressure  occurs.  Note  that  this  happens  before  the  maximum 
eccentricity  is  reached. 

Comparing  Fig.  2(d)  with  Fig.  3(d)  we  see  that  the  peak  pressure  occurs 
at  essentially  the  same  eccentricity  (i.e.,  0.88)  but  that  the  shape  and 
magnitude  of  the  pressures  are  markedly  different.  For  the  PTFE  liner  in 
Fig.  3(d)  the  peak  pressure  is  51  percent  lower  than  the  rigid  liner  of 
Fig.  2(d).  Also  it  is  illustrating  the  effect  of  the  displacement  in  the  high 
pressure  region  producing  a  flattening  shape  along  the  axial  direction. 

Figure  2(e)  represents  the  end  of  the  squeeze  cycle  and  the  beginning  of 
the  separation  cycle.  Now  the  minimum  film  begins  to  increase  which  creates  a 
suction  effect,  an  accompanying  growth  in  the  bubble,  and  a  dissipation  of  the 
pressure  hump  noted  in  Fig.  2(f).  It  is  interesting  to  note  that  the 
eccentricity  is  the  same  (0.88)  for  Figs.  2(d)  and  (f).  The  drastic 
difference  in  the  pressure  plot  is  a  result  of  the  squeeze  velocity  having 
opposite  signs  in  the  two  situations.  The  minimum-film  thickness  continues  to 
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increase  as  it  moves  in  a  clockwise  direction.  The  representation  shown  in 
Fig.  2(a)  is  repeated  as  are  the  representations  for  each  succeeding  orbit. 

Olmensionless  Load  Capacity  Versus  Eccentricity 

Figure  4  represents  the  dimensionless  load  capacity  (w)  versus 
eccentricity  ratio  (e)  for  different  B  parameters.  By  presenting  the 
results  using  dimensionless  load,  (i.e.,  duty  parameter)  the  effect  of 
viscosity  is  normalized  out  and  one  can  focus  on  the  relative  elastic  effects. 
The  predominant  effect  to  note  in  this  figure  is  that  as  B  increases  the  peak 
dimensionless  load  decreases.  These  curves  were  generated  for  two  different 
liner  materials  (i.e.,  steel  and  PTFE).  The  steel  liner  (B  =  0.000025)  is 
represented  with  the  curve  having  the  maximum  peak  load.  The  other  two  curves 
were  generated  with  PTFE  as  the  liner  material.  The  reduced  dimensionless 
peak  load  for  the  PTFE  material  is  a  result  of  larger  deflections  due  to  the 
local  pressures.  Thus  higher  film  thickness  in  the  region  of  higher  pressures 
result  in  a  pressure  relaxation  and  consequent  reduction  in  load  capacity. 

This  is  a  commonly  observed  result. 

CONCLUSIONS 

A  study  of  the  effect  of  bearing  liner  elasticity  on  vapor  cavitation  and 
damper  performance  under  dynamic  conditions  is  presented.  The  motion 
considered  in  this  study  was  that  of  a  journal  undergoing  noncentered  circular 
whirl.  The  bearing  liners  considered  were  steel,  polymethylmethacrylate,  and 
polytetraf 1 uorethyl ene .  The  analysis  was  performed  using  the  Elrod  algorithm 
for  an  isoviscous  lubricant.  A  range  of  viscosities  were  investigated  (i.e., 
0.066  Pa  •  s  to  0.66  Pa  •  s)  to  help  elucidate  the  effect  of  liner 
elasticity.  The  study  revealed  the  following: 

1.  For  the  low  values  of  viscosity  in  this  study  there  was  no  effect  on 
cavitation  due  to  increasing  the  flexibility  by  two  orders  of  magnitude.  As 
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viscosity  was  Increased  the  effect  of  flexibility  became  apparent  by 
Inhibiting  the  Increase  In  bubble  life. 

2.  Increasing  the  viscosity  led  to  an  Increase  In  bubble  life  for  all 
liner  materials  considered. 

3.  As  the  viscosity  Is  Increased,  better  damping  characteristics  can  be 
achieved  using  a  PTFE  liner  over  a  steel  liner  for  a  limited  range  of 
viscosities. 

4.  Increasing  the  flexibility  of  the  bearing  liner  results  in  a  reduction 
and  flattening  of  pressures  along  the  axial  length  of  the  bearing. 

5.  As  the  parameter  B  Increases  the  peak  dimensionless  load  decreases. 
This  parameter  essentially  normalizes  out  the  effect  of  viscosity. 
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TABLE  I.  -  OPERATING  CONDITIONS 


AR ,  m  . 

.  5.0xl0-4 

R,  m . 

L/D  . 

. 1.0 

e  . 

.  0.1  to  0.9 

mj,  rad/s  . 

.  -92.7 

us,  rad/s  . 

. 0 

0,  N/m2  . 

.  1.72xl09 

Pa,  N/m2  . 

.  1.0133xl05 

Pc,  N/m2  . 

. 0.0 

TABLE  II.  -  EFFECT  OF  FLEXIBILITY  ON  CAVITATION 


Li  ner 
material 

Material 

constants 

Vi  scosi ty , 
n(PaS) 

Bubble  life, 
percent  cycle 

Inception, 

ms 

Collapse, 

ms 

a 

E(GPa) 

Steel 

HI 

.30 

206.0 

0.066 

0.025 

22.74 

7.23 

.077 

.029 

■ 

21.74 

9.34 

l^^H 

.099 

.038 

19.44 

13.04 

Ofi  A 

inn 

inn  nn 

1 

.660 

.250 

100.00 

— 

PMMA 

0.40 

2. 

.76 

0.066 

1.160 

76.73 

22.80 

7.03 

PTFE 

0 

.46 

!  0.50 

0.066 

2.960 

76.25 

.077 

3.454 

77.87 

21 .74 

fiBK  [fist 

.099 

4.440 

81.26 

19.64 

6.94 

.264 

11.841 

100.00 

— 

.660 

29.602 

100.00 
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